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Nonclassicality is a key ingredient for quantum enhanced technologies and experiments involving macro¬ 
scopic quantum coherence. Considering various exactly-solvable quantum-oscillator systems, we address the 
role played by the anharmonicity of their potential in the establishment of nonclassical features. Specifically, 
we show that a monotonic relation exists between the the entropic nonlinearity of the considered potentials and 
their ground state nonclassicality, as quantified by the negativity of the Wigner function. In addition, in order to 
clarify the role of squeezing — which is not captured by the negativity of the Wigner function — we focus on 
the Glauber-Sudarshan P-function and address the nonclassicality/nonlinearity relation using the entanglement 
potential. Finally, we consider the case of a generic sixth-order potential confirming the idea that nonlinearity 
is a resource for the generation of nonclassicality and may serve as a guideline for the engineering of quantum 
oscillators. 


I. INTRODUCTION 

At the heart of quantum technologies lies the fact that 
quantum mechanical systems show features, with no classi¬ 
cal counterpart, that may be employed as resources to per¬ 
form specific tasks better or faster than within the classical 
realm HI. In the context of quantum optics, genuine quan¬ 
tum traits of optical systems have led to the emergence of the 
concept of nonclassicality, which characterizes states whose 
effects are not achievable with classical light ia. In particu¬ 
lar, linear models (here intended as systems that induce linear 
transformations of the bosonic mode operators) have attracted 
much attention in the past decades, due to the development 
of experimental platforms able to implement them. In fact, 
the generation of non-classical light, especially in the form of 
squeezed beams, has proven to be an enabling resource for a 
variety of quantum technological applications 0. 

Recently, alternative experimental platforms have been de¬ 
veloped that can also be coherently controlled and described 
as single-mode bosonic systems — including trapped ions i), 
optomechanical systems 0, atoms in optical lattices Ibl, and 
hybrid systems Q ■ The latter naturally embody a playground 
to discuss and test the generation and characterization of gen¬ 
uine quantum features. In particular, they offer the unique 
opportunity to consider nonlinear (or anharmonic) models, 
given that the possibility to host non-linearities is within reach 
of current technologies, in particular for trapped ions 0 and 
optomechanical systems lH. Interestingly, it has been shown 
that the inclusion of nonlinearities in the oscillator poten¬ 
tial uncovers new possibilities to generate nonclassical states 
iHoHni- However, a general framework that encompasses 
these possibilities remains elusive and in particular a thorough 
quantitative assessment of the link between nonlinearity and 
nonclassicality still lacks. 

The aim of this work is to investigate in details the idea that 
nonlinearity is a general resource to generate nonclassicality 
in single-mode bosonic systems constituted of anharmonic os¬ 
cillators. In particular, we will focus on a quantitative assess¬ 


ment of the phenomenon, as we critically consider specific 
quantifiers of nonclassicality and nonlinearity. In fact, identi¬ 
fying proper measures of these quantities is crucial by itself 
and, in particular, different figures of merit exist that cap¬ 
ture different features associated to nonclassicality lfT8l42^ . 
The quantitative connection of the nonlinear behaviour of an 
oscillatory system and the appearance of nonclassicality has 
recently been tested, in the context of nano-mechanical res¬ 
onators, for the Duffing oscillator model ll25l . Here we extend 
such connection and assess its validity for more general sce¬ 
narios, including three families of exactly solvable non-linear 
oscillators and a generic sixth order potential. 

The remainder of this paper is structured as follows. In 
Sec.j^we review the main conceptual tools and establish our 
notation and formalism. First, we introduce and discuss the 
two quantitative measures of nonclassicality that will be used 
throughout the paper, namely entanglement potential and the 
volume of the negative part of the Wigner function. Then, 
we review a recently introduced measure to quantify the non¬ 
linearity of a quantum oscillator 1^ . which in turn is based 
on an entropic measure of non-Gaussianity Eli. In Section 
[n^we analyse the quantitative connection between nonclas¬ 
sicality and nonlinearity for three different nonlinear poten¬ 
tials having an exact solution. We also highlight some differ¬ 
ences between the two measures of nonclassicality (see also 
Ref 1281). In Sec. we address generic (symmmetric) an¬ 
harmonic potential by considering fourth- and sixth-order per¬ 
turbations to the harmonic one. In Sec. lYl we draw our con¬ 
clusions. 


II. NONCLASSICALITY OF A STATE AND 
NONLINEARITY OE A POTENTIAL 

A. Nonclassicality of a Single-mode Bosonic State 

We consider a bosonic system with a single degree of free¬ 
dom, such as a one-dimensional oscillator or a single-mode of 
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a bosonic field. Since we deal with single-mode systems, we 
will not discuss any notion related to entanglement or other 
nonclassical correlations. 

In the most general terms, a quantum state is said to be non- 
classical if the methods of classical statistics fail to describe 
its properties and phenomenology. In the context of quantum 
optics this definition is made precise by using quasiprobability 
distributions in phase space. Here we are not only interested 
in criteria for nonclassicality, but we seek for a quantitative 
characterization. An excellent summary on this topic can be 
found in the introduction of Ref. Il29l . 


where p is the density matrix of the state under scrutiny, |0) 
is the vacuum state at the ancillary port of the beam splitter, 
B is the beam splitter operator, and E[p\ is a suitable mea¬ 
sure of entanglement. Our analysis will be concerned with 
the ground state of a given Hamiltonian model. By dealing 
with pure states, E[p\ can be chosen, with no ambiguity, as 
the the entanglement entropy. This choice corresponds to the 
entropic entanglement potential defined in Ref. 0221 . which 
has been evaluated by truncating the dimension of the Hilbert 
space to a suitable dimensione, ensuring the normalization of 
the state before and after the beam-splitter. 


1. P-nonclassicality and Entanglement Potential 


2. W-nonclassicality 


According to Titulaer and Glauber 0014321 . a quantum state 
is nonclassical when its P function fails to be interpreted as a 
probability distribution in the phase space (see also Refs. 133]- 
l35l ). It has been recently emphasized 0^ that the P func¬ 
tion is the only quasiprobability distribution which can give 
a description that can be completely modelled using classical 
electrodynamics, therefore supporting the idea that to identify 
a classical state it is necessary to use the P function. In this 
paper we dub such fundamental notion as P-nonclassicality. 

The best known way to quantify P-noncIassicality is the 
nonclassical depth E): It quantifies, operationally, the 
amount of thermal noise that is needed in order to render the 
P function of a given state a well-behaved probability dis¬ 
tribution and the corresponding state classical. This measure 
however is not fully suited for our purposes; in fact, while 
we will be interested in establishing a quantitative hierarchy 
of pure non-Gaussian states in terms of their nonclassicality 
at a set nonlinearity of a given potential, it has been proven 
that such states all saturate the nonclassical depth llJTl . i.e. 
they are equally and maximally nonclassical according to this 
measure. 

This obstacle can be overcome by considering the fol¬ 
lowing. It has long been known that coherent states are 
the only pure states that produce uncorrelated outputs when 
mixed by a passive linear-optics device ll38l . Specifically, P- 
nonclassicality has been identified as a necessary condition for 
having entangled states at the output of a beam splitter ll^l40l 
and quantitative relations have been identified between non¬ 
classicality and entanglement ll22ll4TH46l or discord-like cor¬ 
relations more in general ll46l . The idea of quantifying non¬ 
classicality of a single mode state as the two mode entangle¬ 
ment at the output of a linear optic device was introduced by 
Asboth et al. E^ . In particular, it was shown that the opti¬ 
mal entangler is just a (50:50) beam splitter with vacuum as 
an auxiliary state. By restricting to this setup, nonclassicality 
of the input state becomes a necessary and sufficient condition 
for output entanglement. As a consequence, entanglement at 
the output of a beam splitter may be used as a faithful quan¬ 
titative measure of P-nonclassicality. This measure is usually 
referred to as entanglement potential 8 (p) and it is defined as 


8 [ p ]= e \ b { p ^ |0)(0|)Bt 


( 1 ) 


While the P function can be a singular object, the Wigner 
function is always well behaved, even if it can attain negative 
values. Negativities of the Wigner function associated with 
a given state defines the so-called W-nonclassicality, which 
is, however, only a sufficient condition for P-nonclassicality. 
It follows that there are W-classical states which are P- 
nonclassical: displaced squeezed states are a remarkable ex¬ 
ample. The notion of kP-nonclassicality has gained an oper¬ 
ational meaning as follows: the evolution of a system which 
is in a VP-nonclassical state cannot be efficiently simulated 
with classical resources BtI l48l . In order to quantify W- 
nonclassicality we use the volume of the negative part of the 
Wigner function ETTl 

6 = (^Jdxdp\W{x,p)\^ -1, (2) 


where x and p are phase-space variables, and W{x,p) is the 
Wigner function of the state under scrutiny. We will make use 
of the following normalized version of this measure 


V = 


6 

r+5’ 


(3) 


which gives u G [0,1]. 

Let us stress that the W and P-nonclassicality single out 
different quantum features. In particular, the Hudson the¬ 
orem 14^ guarantees that the sole pure states with a posi¬ 
tive Wigner function are Gaussian ones, i.e. squeezed co¬ 
herent states. Hence, there exist pure states that have zero 
VL-nonclassicality (e.g., squeezed states) but non-zero P- 
nonclassicality. In this sense the entanglement potential can 
reveal more detailed features of quantumness, as we will see 
below. Note that measures of VP-nonclassicality based on the 
geometric distance between quantum states have also been in¬ 
troduced ifTSlI^ . 


B. Quantifying the Nonlinearity of a One-dimensional 
Potential Using Its Ground State 

The first idea to quantify the nonlinearity (inteded as the an- 
harmonicity character) of a potential would be defining a dis¬ 
tance between potential functions and the reference harmonic 
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potential. However, this is in general not feasible, since po¬ 
tentials do not need to be integrable functions. A different 
approach follows from the fact that ground states and equi¬ 
librium states of anharmonic potentials are not Gaussian, as 
opposed to those of a quantum harmonic oscillator. We can 
thus choose to quantify nonlinearity by the non-Gaussianity 
of the ground state of a given Hamiltonian model ll26l . The 
measure of non-Gaussianity used for this goal is the entropic 
measure introduced in lIZTl l50l . Here we shall briefly review 
these measures. 


1. Non-Gaussianity of a Quantum State 

The covariance matrix of a single-mode bosonic system 
prepared in a state p is defined as im 

'^jk[p] = {Ri)p{Rk)pi (4) 

where R = is the vector of single-mode quadrature 

operators x and p, and the subscript implies that expectation 
values are calculated over state p. We also dehne the displace¬ 
ment vector X[p] with components Xk[p] = {Rk)p- A Gaus¬ 
sian state has a Gaussian Wigner function. 

To quantify non-Gaussianity of a generic state p, a refer¬ 
ence Gaussian state t should be defined. This is identihed 
as the Gaussian state having the same covariance matrix and 
displacement vector as p. That is 

X[r]=X[p], a[T]=a[p\. (5) 

Non-Gaussianity can now be dehned as the distance be¬ 
tween p and T calculated using, for instance, the quantum rel¬ 
ative entropy 

S'(pIIt) = Tr[p(lnp - Inr)]. ( 6 ) 

We have that S'(p||r) = 0 iff p = r. Although S'(p||r) is 
not symmetric in its arguments, and thus does not embody a 
proper metric, it has been used widely to quantify the distin- 
guishability of two states. This leads to the definition of the 
entropic measure of non-Gaussianity 

4(p) = -S'(pllr) = Tr[plnp] - Tr[plnr] = S{t) - S{p), 

(7) 

where S denotes the von Neumann entropy and, owing the 
way T is defined, we have that S'(r) = — Tr[rlnT] = 
— Tr[plnT . This measure satisfies a series of quite useful 
properties ll27ll : it is additive under the tensor product oper¬ 
ation, and invariant under symplectic transformations, which 
are both very useful for the sake of our analysis. 

The von Neumann entropy of a single mode Gaussian state 
takes the very simple form 

S'(pg) = h(ydet aj , (8) 

where h{x) = (x-l- 5 ) ln(a;-|- ^) — {x— |) ln(a:— \). Thanks 
to this form, the entropic non-Gaussianity becomes 

(5e(p) = h(ydeta^ - S{p), (9) 

which is further simplified for pure states, as S{p) = 0. 


2 . Nonlinearity of a Potential 

We consider a generic potential V(x) and denote with \(j)) 
the ground state of the corresponding Hamiltonian. The hrst 
idea to quantify nonlinearity would be again using the geomet¬ 
rical distance between the ground state of the potential and a 
reference harmonic state, in particular for this purpose the Bu¬ 
res metric has also been employed 1^ . This way of reason¬ 
ing has a downside because we have to choose a value for the 
frequency ui of the reference harmonic oscillator. The most 
natural choice is expanding the potential near its minimum 
and hnding a; as a function of the nonlinear parameters of the 
potential. However, determining this frequency is not always 
straightforward and for some potentials exhbiting more than 
one minimum it may be even be misleading. 

Instead of using a metric we choose to quantify nonlinearity 
using the entropic non-Gaussianity Je, so that the measure of 
nonlinearity is dehned as 

??ng[4^] = <5e (I </>)(</'I) = h(ydetaj, (10) 

this equality holds because the ground state is pure and a is 
the covariance matrix of the ground state (we drop the depen¬ 
dence from the state when obvious). 

This dehnition is more appealing than a geometric one be¬ 
cause it does not require the determination of a reference po¬ 
tential for V(x), but just the reference Gaussian state for the 
ground state of V(x). This makes png independent of the spe- 
cihc features of the potential, since we do not need to know 
the behavior of V {x) near its minimum to compute the refer¬ 
ence frequency. 

Moreover, png inherits the property of the non-Gaussianity 
measure and is invariant under symplectic transforma¬ 
tions EH. This means that (5 ng assigns the same nonlinear¬ 
ity to oscillators which are displaced, rotated in phase space 
or squeezed, which is a reasonable property for a measure of 
nonlinearity. 

III. EXACTLY SOLVABLE NONLINEAR OSCILLATORS 

We now analyze quantitatively the relation between the fig¬ 
ures of merit introduced in the previous Section, considering 
three exactly solvable anharmonic oscillators. 

A. Modified Harmonic Oscillator 

The Modified Harmonic Oscillator (MHO) potential is de¬ 
fined as ll5^ (throughout this manuscript we choose units such 
that h = m = V) 

oP‘xP‘ 

VMHo(a:^) = ^-a/3a;tanh(^a;). (11) 

Here a is a parameter corresponding to the frequency of the 
unmodified harmonic oscillator, while /3 determines the defor¬ 
mation of the harmonic potential. The effects of this parame¬ 
ter on the shape of the potential is appreciated from Fig. 







4 


where 14iHo(a;) is plotted at a set value of a for different 
choices of (3, showing that an increasing deformation param¬ 
eter transforms a harmonic potential into a double-well one 
whose well-depth and separation both increase with j3. 

The normalized wave-function of the ground-state of this 
potential can be found to read ll52l 


0MHo(a^) 


•\/2e 2 cosh(^x) 
yfVl + exp[^Va] 


( 12 ) 


The associated energy is Eq = {a — I3‘^)l2. The covariance 
matrix of such least-energy state can be computed straightfor¬ 
wardly to be 


V{T) 



hNG(r) 


^MHO 


'J_ exp[;3Vc] 

2a ' a^ l+exp[;S^/a] 


|. (13) 


2 l+exp[/3^/a] z 


Its determinant is 


MHO 

det a = - — 


1 T 


2 2T^e^ - -f 1 


4 2 


+ 1)2 


( 14 ) 


with T = ja. Such dependence on r, rather than ot and 
P independently, is common to png = h{det a) and the mea¬ 
sure of nonlinearity based on the Bures distance (for the lat¬ 
ter, we should choose the unmodified harmonic oscillator with 
frequency a as a reference). Both measures of nonlinearity in¬ 
crease monotonically with r. 

The Wigner function associated with ())mho can be written 
in terms of the suitably rescaled phase-space variables q = px 
andp = as 1521 



FIG. 1: (Color online) The MHO potential a = 1 and /3 = 3 (solid 
blue), 2 (dotted yellow), 1 (dashed green) and the harmonic potential 
with unitary frequency and mass (dot-dashed orange). The inset rep¬ 
resent the same graph at a larger scale, where we see the resemblance 
to the harmonic potential. 


FIG. 2: (Color online) Parametric plot of the FF-nonclassicality mea¬ 
sure v{r) versus the degree of nonlinearity for the MHO potential 
and for r € [0.1,6). 


which shows again the key role played by r and, in turn, that 
the non-classicality measure based on the volume of the neg¬ 
ative part of WunoiQjP) determined by such parameter. 

In order to understand how FT-noncIassicality and nonlin¬ 
earity are related to each other, we have studied both quan¬ 
tities against r. In Fig. we report the resulting parametric 
plot, showing that i/ monotonically increases with thus 
supporting the idea that a growing degree of anharmonicity of 
the potential results in increased nonclassicality of the corre¬ 
sponding ground state. 

However, the picture changes significantly as soon as we 
consider P-nonclassicality quantified by the entanglement po¬ 
tential which, as said, can single out more detailed features of 
quantumness. Indeed, at variance with what has been found 
above, such figure of merit turns out to depend on a and P 
independently. The reason for such a difference in behavior 
should be ascribed to the fact that entanglement at the output 
of a beam splitter can be originated either by a non-Gaussian 
input state or by Gaussian single-mode squeezing. In other 
words, nonlinearity is needed to generate W-nonclassicality, 
while P-nonclassicality may be obtained using just squeezing. 

In order to illustrate this clearly, in Fig. we show the en¬ 
tanglement potential and squeezing for the MHO both as a 
function of P for fixed values of r, and as a function of r at 
set values of a. The squeezing in Fig.j^is shown in terms of 
the ratios 


Tx 


,.MHO 

'11 


o MHO 
2(711 ) 


^MHO 

>22 


2(7. 


MHO 

22 


(16) 


with ail = <^22 = 1/2 the variances of position and momen¬ 
tum calculated over the vacuum state of the harmonic poten¬ 
tial. Squeezing is found in the ground state of the MHO for 
either Tx < 1 or < 1. As it is apparent from Fig. the 
behavior of S is rather different from v, and its features may 
be understood looking at squeezing. In particular, we see that 
S grows when the ground state exhibits squeezing. 
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(a) 


(b) 



FIG. 3: (Color online) Entanglement potential and squeezing for 
the MHO. In panels (a) and (b) [(c) and (d)] we plot (dashed 
blue curve), Vp (dotted orange curve), the entanglement potential (P- 
nonclassicality) £ (black dots), and FF-nonclassicality zz (red dots) 
against f5 [r] for r = 1 and r = 3 [a = 3 and a = 5]. Squeezing is 
observed for either < 1 or rp < 1 (i.e. variances of the perturbed 
ground state below the values of the vacuum state of a harmonic os¬ 
cillator). 


B. Morse potential 

The Morse potential has been introduced as an approxima¬ 
tion to the potential energy of diatomic molecules as it pro¬ 
vides a better description of the vibrational structure than the 
(quantum) harmonic oscillator ll5^ . The form of the potential 
is 

= (17) 

where x is the distance from the minimum of the potential, the 
parameter D > 0 determines the depth of the well, while a 
controls its width. Expanding the two exponentials for a —0 
at fixed D we get the harmonic limit, which is an oscillator 
with frequency cum = \/2Da. The potential is plotted in 
Fig.|^for different values of the parameters. 

The Schrodinger equation associated with this potential 
can be solved analytically, the energy eigenvectors being la¬ 
belled by two quantum numbers, which we label here N and 
V. The first is r elate d to the parameters of the potential as 
N = -1/2 -F y/Wja. The second, which can take values 
V = 0,1,2,.., [7VJ, counts the number of anharmonic excita¬ 
tions of the system. As we want at least one bound state, we 
require > 0. We thus have the constraint a < 2\j2D. The 
limiting case where we have just one bound state (the ground 
state) is achieved for a 2\/2D. The wave-function of the 
ground state is 

Mx) = {2N + (18) 

with associated energy E = —aN‘^/2. The behavior of the 
nonlinearity of the Morse potential can be understood by look¬ 


'/m(x) V„(x) 



FIG. 4: (Color online) The Morse potential Vm{x) for D = 1 and 
a = 1 (solid blue), 2 (dashed orange), 3 (dotted green). The inset 
shows the potential for a = 1 and D = 1 (solid blue), 2 (dashed 
orange), 3 (dotted green). 


ing at the form of the potential in Fig. as opposed to the 
harmonic one ll^ : For any fixed value of D [a] we expect an 
increase [decrease] of nonlinearity for increasing a {D}. 

The covariance matrix associated with the ground state in 
Eq. ([^ is 






(19) 


where is the polygamma function = 

logr(z), T(Z) being the Euler Gamma function. The 
determinant of this correlation matrix, and thus the Bures dis¬ 
tance from the reference harmonic oscillator, depend just on 
N or, equivalently, on the combination \J2E ja. In this case 
both measures of nonlinearity are monotonically decreasing 
functions of N. 

The Wigner function for the ground state of the Morse po¬ 
tential reads as follows ll54ll 

W ^\ - „-2NaXT^ froAT , T\„-ax\ 

^m\XtP) 7rr(2A^) ^ 2ip/ct (( 2 A + Ije j, 

( 20 ) 

where (z) is the Macdonald function of (non-integer) order 
7. In order to calculate the measure of nonclassicality iz, we 
rescale the phase-space variables to q = ax and P = g, and 
evaluate 


JJ da;dp |lTM(a;,p)| = 

2g-2tVq 


dp dp 


TrT(2N) 


{2N + 1)^^ K_2^p{i2N + l)e-’>) 


( 21 ) 


which shows that the only relevant parameter is N. The nu¬ 
merical integration of Eq. © is challenging and was carried 
out with the aid of the CUBA libraries ll55l . The degree of FU- 
nonclassicality iz is found to monotonically decreases with N, 
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and the parametric plot of nonclassicality versus nonlinearity 
in Fig. [^reveals a monotonic behavior, strengthening the link 
between such features and reinforcing the idea that nonlinear¬ 
ity might play the role of a catalyst for nonclassicality. 


v(N) 
0.4 ■ 



FIG. 5: (Color online) Parametric plot of the FF-nonclassicality v 
versus the degree of nonlinearity rj for a Morse potential with D = 1 
and a £ [0.15, 2.7], i.e. N € [0.0238, 8.928[. 


The situation regarding the entanglement potential is com¬ 
pletely analogous to what we found for the MHO, as it de¬ 
pends on both parameters. In Fig. |^we report the same kind 
of graphs, both with N fixed and a hxed, which show that the 
behavior of £ is explained by the squeezing of the state. 


(a) (b) 






FIG. 6: (Color online) Entanglement potential and squeezing for 
the Morse oscillator. In panels (a) and (b) [(c) and (d)] we plot 
(dashed blue curve), Vp (dotted orange curve), the entanglement po¬ 
tential (P-nonclassicality) £ (black dots), and FF-nonclassicality v 
(red dots) against a [l/A^] for N = 1 and = 5 [a = 1 and 
a = 3]. Squeezing is observed for either < 1 or rp < 1 (i.e. vari¬ 
ances of the perturbed ground state below the values of the vacuum 
state of a harmonic oscillator). 



FIG. 7: (Color online) The Posh-Teller potential with Apt = 1 and 
a = 1/2 (solid blue), 1 (dotted orange) and 3 (dashed green). 


C. Poscbl-Teller potential 

The modified Poschl-Teller potential (PT) is defined as 

VpT(a;) = —ApTCOsh“^(aa:), (22) 

where Apt > 0 is the depth of the potential and a is con¬ 
nected to its range. The harmonic limit is obtained at hxed 
Apt for a —> 0 and the frequency of the reference har¬ 
monic oscillator is wpT = y/2A^a. As for the Morse 
potential, we have a quantum number s that labels the en¬ 
ergy eigenstates and counts the anharmonic excitations. It 
is related to the parameters of the potential through the re¬ 
lation Apt = |a^s(s + 1). Therefore, the request for 
the existence of at least one bound statefranslates into s = 
I ^—1 -F y/l -F SApt/q!^^ > 0. Fig. m shows the depen¬ 
dence of the PT potential on the position coordinate. 

The ground state of the system reads 


(j^vrix) = \ ^^ cosh ^(ax), (23) 

7r4 y r(s) 

with associated energy E = 

Differently from the previous cases, the covariance matrix 
of the ground state is rather involved and will not be reported 
here. In line with the case of the previous two anharmonic 
potentials studied here, its determinant depends only on s (or, 
equivalently, on App/a^). Again, both png and the Bures 
nonlinearity are monotonically decreasing function of s only. 

The Wigner function of state ^pT(a:) in Eq. ( |2^ is known 
analytically for the case of Apt = ll52]| . In this case, the 

measure v is an s-dependent constant, as it can be seen by 
rescaling the relevant variables as p' = ^, x' = ax, y' = ay 
and evaluating the integral 

Wvy{,x,p) = JAy (t)^{x- |)</>pT(a;-F (24) 

which embodies the definition of Wigner function. 
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(a) 


(b) 





FIG. 8: (Color online) Entanglement potential and squeezing for the 
PT oscillator. Panels (a) and (b) [(c) and (d)] show (dashed blue 
curve), Tp (dotted orange curves), and the entanglement potential £ 
(black dots)against a [1/s] for s = 1 and s = 5 [a = 1 and a = 3]. 


As for the entanglement potential, this turns out to depend 
on both a and s. Plots similar to those valid for the MHO 
and Morse potential are presented in Fig. (without the W- 
nonclassicality tz). 
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FIG. 9: (Color online) Contour plot of the overlap between the per¬ 
turbative ground states of Eq. {16\ and the numerically calculated 
one (for ui = 1). 


tion theory ll5^ . The state takes the form 


IV. OSCILLATORS WITH POLYNOMIAL 
PERTURBATIONS 


3 


l'0) = X!'^2n|2n), 

n—0 


(26) 


So far we have studied exactly solvable potentials with two 
parameters and revealed a common behavior: the nonlinearity 
and the FF-nonclassicality i/ have the same behavior and de¬ 
pend just on a single effective parameter. On the other hand, 
the entanglement potential carries a dependence on both the 
parameters and its different behavior may be understood in 
terms of the squeezing of the state. 

Now we want to address the case of a generic two- 
parameter perturbation, so we study a physical system com¬ 
posed of a one-dimensional harmonic oscillator with perturba¬ 
tions proportional to and x® respectively. The Hamiltonian 
of this system thus reads 

H = uP'x^) -\- £4X^ + cgx® . (25) 

As the model is not exactly solvable, the properties of the sys¬ 
tem will be studied using perturbation theory. We notice that 
Eq. ( |25l l may also serve as an approximation for any symmet¬ 
ric (even) potential. We do not consider odd powers of x, as 
they give rise to known problems in the convergence of the 
perturbative series. A remark is in order: terms proportional 
to X and to x^ could in principle be treated in a perturbative 
way as well. However, they do not give rise to truly anhar- 
monic behavior, and will not be considered in this context. 

In order to get insight into the ground states for these 
Hamiltonians we use hrst-order time-independent perturba- 


where \k) denotes a Fock number state of the harmonic oscil¬ 
lator and the coefficients 7 ^ are in given Appendix ]A 1| 
First-order perturbation theory returns the ground as a h- 
nite superposition of Fock states, which makes the Wigner 
function and the nonlinearity easy to compute. In order to 
assess the validity of the first order approximation, we com¬ 
pare such ground state to the state obtained by numerically di¬ 
agonalizing the Hamiltonian of the system within a truncated 
Fock space of suitable size. Convergence of the results of 
such numerical calculations appear to be ensured by using 61 
harmonic levels. The corresponding ground state \(j>) is then 
compared to \tj}) using the state hdelity \{(j)\ijj)\'^. In Fig.we 
present a contour plot of the overlap as a function of both €4 
and eg. For values of £4 up to 0.1 and £5 up to 0.03 the fidelity 
is at least « 0.976. 


A. Nonclassicality and Nonlinearity 

From the perturbed ground state in Eq. ( |2^ we compute 
the nonlinearity of the perturbing potential. The covariance 
matrix associated with \ip) can be thus written as 

(l-r2(a")+2(ata)-4(a)^) „ \ 

0 f(l + 2(ata)-2(a2))y’ 

(27) 
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v{£4.£e) 



FIG. 10: (Color online) Random scatter plot of the W- 
nonclassicality 1 / versus the nonlinearity ? 7 ng for the perturbed har¬ 
monic oscillator when both parameters £4 and eg are are varied in the 
range given in Fig.[^ 1000 random points were generated. The dark 
blue curve below the points represents eg = 0 , while the one above 
the points is the curve for eg = 0.03. The light red curve below the 
points is the one for £4 = 0 . 1 , while the one above the points is for 
£4 = 0. 


£(^4.^6) 



FIG. 11: (Color online) Random scatter plot of the entanglement 
potential (P-nonclassicality) £ versus the nonlinearity r^NO for the 
perturbed harmonic oscillator when both parameters £4 and £g are 
varied in the range given in Fig.[^ 10000 random points were gener¬ 
ated. The dark blue curve above the points represents £g = 0, while 
the one below the points is the curve for £g = 0.03. The light red 
curve above the points is the one for £4 = 0 . 1 , while the one below 
the points is for £4 = 0 . 


with a and the annihilation and creation operators of the 
oscillator and 

(a'l') = (a) = 0, 

(a’I'a) = 2|72p +4|74p -F6|74p, 

(a^) = V27270* + 2V37472* + V^7674*, 

(at2) = (^2)* = V2jol2* + 2^37274* + 'MiaIq* ■ 

(28) 


An explicit calculation shows th at the determination of 17 ^°*, 
and in turn the nonlinearity /i(-\/det a), depends on both the 
perturbative parameters and on the frequency w. No single¬ 
parameter rescaling can be identihed in this case, thus entail¬ 
ing the double-dependence highlighted above, which is passed 
to the FF-nonclassicality v [cf. Appendix ] A 1| . 

As our aim is to highlight the role played by the perturbative 
parameters, in the remainder of our analysis we set a; = 1 and 
generate random pairs of values (£ 4 , eg) (within the appropri¬ 
ate range of validity of the hrst-order perturbative approach 
discussed above) that are then used to compute both the non- 
classicality and nonlinearity indicators. 

The results shown in Figs. and 11 showcase a non¬ 
monotonic relation between nonlinearity and nonclassicality; 
the points corresponding to the randomly taken pairs of val¬ 
ues for the parameters are distributed within a (narrow) region 
comprised within four curves, each associated with an ex¬ 
tremal value of £ 4 , 6 . Nonclassicality and nonlinearity are thus 
strongly dependent on the details of the system under consid¬ 
eration and are, strictly speaking, non equivalent notions. On 
the other hand, the regions in Fig. and 11 are concentraded 
enough to suggest that the intuitive link between such two fea¬ 
tures is, in fact, correct. Notice also that upon fixing the value 
of one of the parameters (either £4 or £ 5 ) the behaviour of both 
nonclassicality measures becomes monotonic with nonlinear¬ 
ity. 


Digging into the details of the phenomenological behavior 
identihed by our analysis, it appears that FF-nonclassicality 
(Fig. 10 1 is favoured by the -like nonlinearity. On the other 
hand, P-nonclassicality, appears to beneht from a i^-type of 
nonlinear effects: in Fig.[^the roles of the dark blue and light 
red curves are inverted with respect to Fig. showing that, 
after choosing the parameters £4 and £g in such a way that the 
entropic nonlinearity is hxed, the ground state obtained with 
the maximum value of £4 generates more entanglement than 
any other one. 


V. CONCLUSIONS 

We have addressed in details the role played by the nonlin¬ 
earity of anharmonic potentials in the generation of nonclas¬ 
sicality in their ground states. In particular, we have shown 
that nonlinearity plays a crucial role in the generation of W- 
nonclassicality, while P-nonclassicality may be also obtained 
by potential inducing just squeezing. 

Our results support the expectation, put forward in 
Ref ina, that the nonlinearity of a potential is quantitatively 
related to the nonclassicality of its ground state and thus the 
former feature may be regarded as a resource to generate the 
latter one. The strict validity of such expectation, which ap¬ 
pears to be conceptually quite intuitive, is however strongly 
linked to the specihc details of the Hamiltonian model being 
addressed. Anharmonic potentials that can be reduced to a 
single-parameter dependence give rise, in fact, to a monotonic 
relation between nonlinearity and W-nonclassicality. Such a 
correspondence breaks down for effectively multi-parameter 
potentials: set values of nonlinearity bound the possible de¬ 
grees of nonclassicality of the ground state of a given anhar¬ 
monic potential, albeit without determining it unambiguously. 

Our investigation opens up a series of questions, all linked 
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to the effective role that non-harmonic oscillators might have 
in the quantum technology arena; it would be interesting, for 
instance, to investigate whether the enhanced non-classicality 
achieved, in general, for a non-null degree of nonlinearity is 
accompanied by an equally enhanced degree of coherence in 
the ground state of the oscillator. Equally interesting is the 
question on the actual use that can be made of the sought non¬ 
linearity in protocols of practical quantum estimation. 
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Appendix A: Calculations for the harmonic oscillator with perturbations 
1. Perturbative states 

The matrix elements of the two perturbations on the basis of the energy eigenstates of the unperturbed system, which in this 
case are the Pock states \n), are 


(n|a;^|n) = 
{n\x‘^\n -I- 4) = 
{n\x'^\n + 2) = 


Q-n? -I- 6n -I- 3 

~ 4w2 ' 

\J{n+ l)(n -f 2)(n -f 3)(n -I- 4) 

^ ’ 

(4n + 6) yjn{n — 1) 

4^ ’ 


for x^ and 


{n\x^\n + 6) 
{n\x^\n + 4) 
{n\x^\n + 2) 


5(4n^ -I- -f 8n -I- 3) 


8w3 


{n + l){n + 2){n -f 3)(n -f 4){n -f 5)(n -f 6) 


8a;3 


3(2n -f b)yj{n + l)(n -f 2)(n -f 3)(n -f 4) 


8w3 


15 (n^ -f Sn -f 3) 1/(71 -f l)(n -f 2) 


8a;3 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 


for and all the other elements are zero apart from the symmetrical ones (i.e. {n\x‘^\n + k) = (77-|-fc|x^|n) and {n\x^\n + k) = 
(n -I- k\x^\n)). 

The formula for the perturbed ground state is the following 




(A8) 


where Vnk = and V stands for the perturbation e 4 X^ + cqX^. Using this formula and the matrix elements found in 

the last section we readily find the normalized ground state ((26]l, its coefficients are the following 


1 

C 

74 = -70 


7o f 45e6 ^ 

^/2 V 4a;3 


/3/15e6 ^ e4\ 

V 2 2 a ;3 ) 


(A9) 


76 


= - 1 / 57066 , 
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where the normalization constant C is 


C = 


V'a;2(96a;6 + llTe^) + 945a;e4e6 + 2055i| 


4\/6w^ 


(AlO) 


2. Wigner function 


It is convenient to express the Wigner function as W{a) = ^Tr ll57l l58l where D is the displacement 

operator D{^) = exp(^a^ — ^*a). The expectation values of D on Fock states are given by 


. if n'>n 

{n'\Diz)\n) = { ^ 

if n>n' 


(All) 


where L^n\x) are the associated Laguerre polynomials. The Wigner function then becomes 


W{z) = [7o2Lo(4|zn + ^2^L2{4\z\'') + 

TT 

+ 76^A6(4|zn + 4^27072 Re{z‘^)Ll{‘i\z\'^) 

1 fi *^2 

+ ^7074Re(z^)T^(4|zn + 

+ ;^7274R.e(z^)L4(4|z|^) + ^^7275 Re(z'‘)L^(4|z|2) 


8 


730 


7476Re(7)L74|^n], 


(A12) 


where the coefficients are given by 
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